One-Dimensional XXZ Model for Particles Obeying Fractional Statistics 
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We define one-dimensional particles as non-abelian representations of the symmetric group Sn- 
The exact solution of an XXZ type Hamiltonian built up with such particles is achieved using the 
coordinate Bethe Ansatz. The Bethe equations show that fractional statistics, effectively, accounts 
for coupling an external gauge field to an integer statistics' system. 



PACS Numbers: 71.10.Pm, 71.27.+a, 75.10.Jm 

Physical behaviour of quantum systems is deeply af- 
fected by the statistics of the constituting effective de- 
grees of freedom. Quasi-particles and quasi-holes in con- 
densed matter physics may obey statistics interpolating 
between fermionic and bosonic behaviour. Examples are 
the excitations of two-dimensional electron systems ex- 
hibiting Fractional Quantum Hall effect These exci- 
tations are called anyons. They have been a subject of in- 
tense study also in connection with superconductivity ||] 
and superfluidity j| . Fractional statistics of such parti- 
cles arises from the trajectory-dependence of the particle 
exchange procedure in the two-dimensional configuration 
space. This feature makes the concept of anyons purely 
two-dimensional. The Fock space formulation of anyon 
operator algebras takes into account these characteristics. 
The creation and annihilation operators (introduced as 
Jordan- Wigner transforms of usual fermions on a two- 
dimensional lattice |Q or as unitary representations of 
the diffeomorphism group of R 2 || ) obey deformed com- 
mutation relations if the exchange involves anyons at dif- 
ferent spatial positions (see Appendix). 7V-anyon-states 
are abelian representations of the braid group Bn ||] 
(whereas bosons and fermions furnish, respectively, the 
identical and alternating abelian representations of the 
symmetric group Sn)- These features make anyons dif- 
ferent from g-oscillators, the latter providing a realization 
of Gel'fand-Farlie quantum group, which is a local defor- 
mation of the Weyl-Heisenberg (bosons) or Clifford alge- 
bra (fermions) 0. The path dependence implies that the 
one-particle state is inextricably related with the com- 
plete state of the many body configuration. This intrin- 
sic non-locality makes anyon physics very difficult. Even 
statistical properties of a free anyon gas are only partially 
established using the virial expansion || . 

Haldane || formulated the notion of fractional statis- 
tics without any reference to the spatial dimension D. 
The generalized Pauli principle is expressed in terms of 
the reduction of the single-particle Hilbert space when 
particles are added to a many body system keeping 
boundary conditions fixed. Another way to introduce 
dimensionality-independent fractional statistics has been 
formulated in Ref. JlO[ where quons have been intro- 
duced. Quons' fractional statistics results from the 
"superposition" of statistical properties of bosons and 



fermions pi]] . In D > 2, this is consistent with spin- 
statistics theorem. In 2D, Haldane particles and quons 
capture the essential features of anyons (||| ■ 

Recently, an outgrowing interest has been devoted to 
generalized statistics in one dimension. A specific way 
to introduce D = 1 fractional statistics has been pro- 
posed in connection to the quantization of the solutions 
of the Calogero model There, the potential 1/x 2 

is interpreted as "statistics interaction". The same no- 
tion of fractional statistics applies also to anyons in a 
strong magnetic field that restricts the allowed energies 
to the lowest Landau level. The anyon gas, then, is de- 
scribed by an effective field theory on a ring where the 
dynamics of particles is one-dimensional [fl3[ . It is worth- 
while noting that such one-dimensional particles obey 
fractional statistics, but they are not "true anyons" since 
in D =/= 2 trajectories in the particle configuration space 
have no meaningful braiding property. Instead, nonlo- 
cal "deformations" of the commutation relations furnish 
non-abelian representations of the symmetric group Sn- 

In this paper, we deal with particles in D = 1 that 
preserve the intrinsic non-locality of two-dimensional 
anyons, but which are still representations of Sn- This 
representation is no longer abelian. The second quan- 
tized formalism and the Fock space representation is de- 
veloped. The XXZ model for such particles is formu- 
lated and solved exactly using coordinate Bethe Ansatz 
(BA) Q in D = 1. 

For D = 1 we define a set of creation/annihilation op- 
erators {// , fi} for a spinless particle at site i. They 
obey the deformed relations 



fjfk + Qj.k fkfj = S Jtk 

fjfk + qj'jfkfj = o 



(i) 

(2) 



where q~l = (<Zj,fc) _1 . Since the operators are path- 
independent in D 2 (compare Appendix), (Q) and (|^) 
have to constitute a representation of Sn, and not of Bn- 
This is ensured by the "consistency relations" 



9i,k = ?*, 



(3) 



(4) 
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Such a representation is non-abelian ra. For j — k, (g) 
gives qj t j — ±1. Hence (Eh, (|^) are an extension of anyon 
commutation relations pjjq] to D ^ 2. In contrast to 
the true anyonic case (see Appendix), Qj tk has no rela- 
tion with the configuration space geometry, but is a free 
"external" parameter. 

Relations (|l|) and (|^) are formally analog to quon com- 
mutation rules [jllj . Note that the deformation parame- 
ter here depends on two indices (j,k), whereas it does 
not in quon commutation rules. Without this index- 
dependence, relation (||) directly implies q 2 — 1. As 
a consequence, quons obey integer statistics in ID if q 
is a C-number [fi"l| (as an operator, it has eigenvalues 
±1). If q hk = ±11 V (j,k), then (|), (|) describe spin- 
less fermions or bosons, respectively. However, for appli- 
cation we choose qj t k being C- numbers and qjj — 1, see 
@, which implies the Pauli exclusion principle, as for 
spinlcss electrons or hard core bosons. 

Relations (||), (||) hold if qj tk is an operator commuting 
or anticommuting with both, fj and ft. For this reason 



we add this as a postulate |T(| 



Ifl 



Qj,k] = [fj , Qj,k] = 



(5) 



The introduction of two indices for the deformation pa- 
rameter allows the construction of consistent commuta- 
tion relations even for qj.k being C-numbers. We make 
use of this possibility in 

To develop a Fock representation of the algebra ([!]), 
(||), we take Vj = fjfj as number operators. Relations 

©i (§)' an d @ yield commutators of Vj and fj, fj being 
unaffected by the deformation parameter q^k'- 

[vj,v k ] = , [vj,fl\ = 5j. k fl , [vj,fk] = Sj,kfk ■ (6) 

Moreover, the property qjj = 11 implies that number op- 
erators are idempotent: (vj) 2 = fj. Because of (||) the 
one-particle Fock representation of the algebra ([!]) , (||) is 
unaffected by Qj k- Instead, the action of /, , vi on the 
TV-particle state \ni . . . n^) is deformed according to 



fl\m..n N ) = 



fl\nt..n N ) = 



i-i 

n*. 

,/c=i 

7-1 

n 

k=l 



\nx..ni - l..n N ) , 



Qk,l 



\n\..ni + l..n N ) , 



vi\n 1 ..n N ) = ni \nx..n N ) , 



(7) 



where ni 6 {0,1}. Equations ([?]) generalize the cor- 
responding relations fulfilled by integer statistics par- 
ticles |l7j characterized by Ylk=i %k = (i) i_1 (f° r 
fcrmions/bosons). An explicit realization of the oper- 
ators fj in terms of spinless fermionic operators ctj is 



fj 



exp (— i $;n;), where $z are hermitean oper- 



ators commuting with fermionic degrees of freedom. By 
direct calculation, the relations (Q), (||) are obtained by 



setting qjk = exp [i($k — $j)]. This realization has been 
suggested in the Ref. |ll| where = Pk, Pk being mo- 
menta of a phononic bath coupled to fermionic degrees 
of freedom. 

In the following we consider the ID anisotropic Heisen- 
berg model (XX Z model) of spinless fermions 



The /-operators obey relations ([j]), (g) with qj.k defined 
in close analogy with anyonic relations (see Appendix) 



qj.k 



j > k 
j = k 
j <k 



(9) 



We point out that the postulate (Q) is fulfilled since qj >k is 
a C-number for arbitrary, fixed (j, k). Relations (||) and 
(^) imply that q is on the unit circle, i.e. \q\ = 1. Addi- 
tionally, Periodic Boundary Conditions (PBC) fi+L = fi 
are chosen, where L denotes the period. The parameters 
t and U are the hopping amplitude and the Coulomb 
interaction strength, respectively. 

It is worthwhile mentioning that (||) implies fixing an 
order on the infinite periodic chain. This order can 
only be defined locally on the manifold S 1 , which needs 
two charts for its description. We choose two "charts" 

Ci = {h,-,3N} and C 2 = {j/v-t, -JnJi, ••, Jn-I-i} 
on the ring thought as a discrete subset of S 1 . On each 
chart, the given order is well defined by interpreting them 
as ordered sets. The intersections between C\ and C2 are 
Un-i ■ •• jiv} and {j x . . . jjv-i-i}- In such sets, the or- 
ders defined on C\ , Ci are identical. Now, in D = 1 , only 
nearest neighbour (n.n.) exchanges can take place. Thus 
qj k .j, is connected to the n.n. exchange jk ji- On the 
chart Ci, where j± < jjy, the exchange jx <-> j\ is not 
a n.n. exchange. To allow for n.n. hopping <-> ji, we 
must use C'2 on which jn < ji- This implies qj N ,j t = q 1 - 
The picture depicted above is equivalent to fixing a pe- 
riod Pq = { ji , ... , ji+L } on the infinite periodic chain. 
Consistency of the PBC with this induced order is given 
if the results are independent of Pq. In the sequel it will 
be seen that this condition is fulfilled. 

The correspondence between and the deformed 
anisotropic Heisenberg model can be established by 



/j, Sf> = fj and S) z) = 1/2 



' 3 5 3 

c(±) 



On-site, the 



operators Sj ,Sj generate the fundamental represen- 
tation (spin s — 1/2) of sit (2), but for j ^ k 



qi + ) q(-) 

Dj , O k 

q(-) <?(") 



= (l + qj , k )S^Sp -Sj, k , 



s 



{z) Q (±) 



{i + qj.k)s\ 
. 



^SP- hk 



(10) 
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We now show that the XXZ model is exactly solvable 
by means of the coordinate BA. The general iV-particle 
state on a chain with L sites can be written as 



I*) 



E 

i<ji <•■■<:/«<£ 



1pi3uh,:;jN) fl fj 2 ■ ■ ■ fj N 1 0} 



(11) 



The action of Hxxz on \ty) , i.e. the eigenvalue equa- 
tion, then reads 



N-l 



1=1 



,3n) 



N 



i=i 



= Eip(ji,...,j N ) 



where = Yl^Ji 1h ,j s 1j s j&i- Using @ as definition 
for qj,k, it follows that 



(12) 



QV = Q\ +) = i • 



(13) 



Hence, exactly the same eigenvalue equation as in the or- 
dinary XXZ model is obtained. Consequently, inserting 
the BA 



V'O'i, -3n) 



( N 

^ A(w) exp \i V 



rn=l 



(14) 



the energy E and the two-body scattering matrix 
S(k, k!) = — exp [— i9(k , k')] for the "deformed" model 
are identical with the known terms occuring in the usual 
XXZ model |lq|. Imposing PBC, however, yields 



,3N,3l) 



' N 
.1=2 

{-q) N -H{3u32, 



ip(juh,---,3N) 



,3n) 



(15) 



Equation (jT||) shows that the fractional statistics pro- 
duces a twist in the PBC that modifies the periodicity of 
the Bethe wave function. 

Since the twisting factor q N_1 does not depend on ji, 
the starting point of the chosen period Po, the bound- 
ary condition is consistent with our choice of <jj t k. The 
twist g^ -1 does not depend on the particles' configura- 
tion, but on the number of particles only. This is crucial 
because configuration-dependence of the twisting factor 
would destroy the solvability of the model, since it modi- 
fies the structure of the exponential functions in imposing 
PBC on the BA wave function (making it impossible to 
extract from ( p"5| ) a relation for the amplitudes A(n)). So, 
the coordinate BA solvability of the model (with PBC) 
demands a careful choice of the functional form of q^h- 
The choice = exp [iS (j — k)] for the XXZ model, 



for instance, leads to the same structure of the S matrix, 
but produces incompatible boundary conditions. 

Since \q\ = 1, q = exp [i arg(q)]. So a phase shift by 
multiples of arg(g) occurs in the BA wave function on 
the r.h.s. of eq. (|l5|). The Bethe Equations (BE) are 
obtained as 



arg([-g 



N-l\ 



27r/, 



N 

E 

rn — 1 



0(k 3 , k m ) , (16) 



where Ij G Z. In the fcrmionic case, one obtains kjL 



N-l' 



Eto=i ^(^j ' k m ): whereas the hard core 



bosonic case yields kjL = 2nlj — J2 m =i ^(^j » ^m)< I n au 
cases, Ij G Z. Equation ( |ilf ) differs from the BE for the 
ordinary XXZ model in the additive term arg([— g]^ 1 ). 
This term has its origin in the fractional statistics of the 
particles, and vanishes for integer statistics. Equations 
dl6| ) were obtained for the ID XXZ model on a ring 
threaded by an external magnetic flux |n],^(J. We re- 
cover the BE of identifying $ = a = argQ-g]^ -1 ), 
$ being the magnetic flux in units of Tic/e. 

The limit fcry —> tt ~ fi, where cos(^t) = U/2t for 
\U\ < 2\t\, in the 7V-th equation of ( pf ) relates the statis- 
tics with the ratio U/2t for the ground state, 



7T (L - N - 2I N + 1) - n (L - 2N + 2) 



(17) 



At half filling (N/L = 1/2, I N = (N - l)/2) the energy 
and the total momentum, P = ^2iLi ki/L, of the ground 
state are affected by the statistics' factor a pp| , p^[ as 
follows: 



E (a)-E {0) 



it sin(/x) 

4/x(7T — (jL)L 



P = - a 
2 



(18) 



where Eq(0) denotes the ground state energy of the un- 
deformed XXZ model. The same structure of BE has 
been obtained in |2l],|22| . In |22| , the two species of par- 
ticles, up-spin (er = +) and down-spin (a = — ), have 
a dynamics governed by two distinct XXZ^ Hamil- 
tonians coupled only via a local gauge field, included 



{<*) 



in XX Z^ by a Peierls-like substitution t ^ tW, 

where Wm = exp [icr^2f = i ot-m-ini,-v\ (determined by 
the position of all particles of opposite species); ai G HI, 
ctm+L = ot m . A comparison of the BE in p2] with jl^ ) 
shows that our deformation parameter q can be inter- 
preted as the "global" coupling constant of the gauge 
potential by setting a ~ X)m=i a " 1 ~ i - Vice versa, such 
an interaction produces statistics transmutation. In this 
sense, our deformed XXZ model belongs to the same 
class of integrable models introduced in p2| . 

In conclusion, we have given a formulation of frac- 
tional statistics in one dimension realized by an anyonic- 
type deformation of the second quantized commutation 
rules. Coordinate BA solvability of the deformed XXZ 
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model demands a proper choice of the functional form 
of 1j,k- The statistics we have chosen in the present pa- 
per preserves the Yang-Baxter equation as well as the 
BA solvability of the undeformed model. The resulting 
BE are, however, modified. They show that fractional 
statistics plays the same role as a gauge field coupled 
to the undeformed model. Systematic investigations of 
fractional statistics seem interesting for at least two rea- 
sons. First, fractional statistics may be an alternative ap- 
proach to handle complicated interactions between parti- 
cles obeying integer statistics. Such interactions could be 
modelled deforming the particles' statistics. Second, the 
study of "compatible statistics" could be relevant in order 
to find integrable Hamiltonians characterized by braided 
Yang-Baxter Equations (YBE) Such a feature of 

the YBE could be closely related to actual braiding of 
particles in two dimensions. 

A further development of the present approach is to 
take spin into account. A preliminary analysis of the 
"deformed" Hubbard model |L5| shows that fractional 
statistics modifies the S matrix; the R matrix obeys a 
braided YBE. We will report on this subject in a forth- 
coming paper. 
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APPENDIX A: 

Here we summarize the commutation properties of two 
dimensional anyons The creation/annihilation op- 

erators obey 

6 f (x c ) b(yc) + <z(x c , yc)b(yc) b*(xc) = 5 
H x c) b(yc) + 9 _1 (xc,yc) b(y c ) b(x c 

The operators tf(xc) (6(xc)) create (annihilate) an 
anyon at site xc = {x\,X2)- C denotes the path run- 
ning from +oo to xc keeping x 2 constant. The rela- 
tions above hold if xp > yc; in the case xp < yc, 
they are satisfied substituting q «-> q^ 1 . Note that 
xc > yc <=> {x2 > Vi V xi > ?/i (if x 2 = y 2 )}- The 
function q(x.c,yc) = <z(l x c —yc\) can be simplified (see 
e.g. Q) to q — e nv (y Ed), where v denotes the statis- 
tics. If two anyons are at the same position xc = yc. 



then 3 = 1. Otherwise, the standard bosonic or fermionic 
algebras are deformed by the parameter q. 



xcyc ) 
. (Al) 
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